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We predict and numerically observe strong periodic orbit 
effects in the properties of open quantum systems with a 
chaotic classical limit. Antiscars lead to a large number of 
exponentially narrow resonances when the opening is located 
on a short unstable orbit of the closed system; the probability 
to remain in the system at long times is thus exponentially 
enhanced over the random matrix theory prediction. The dis- 
tribution of resonance widths and the probability to remain 
are quantitatively given in terms of only the stability matrix 
of the orbit on which the opening is placed. The long-time 
remaining probability density is non-trivially distributed over 
the available phase space; it can be enhanced or suppressed 
near orbits other than the one on which the lead is located, 
depending on the periods and classical actions of these other 
orbits. These effects of the short periodic orbits on quantum 
decay rates have no classical counterpart, and first appear 
on times scales much larger than the Heisenberg time of the 
system. All the predictions are quantitatively compared with 
numerical data. 



I. INTRODUCTION 

Wavefunction scarring, the enhancement or suppres- 
sion of quantum eigenstate intensity along an unstable 
orbit of the corresponding classical system, is a fascinat- 
ing and generic property of quantum chaotic behavior. 
Along with dynamical localization, it is one of the strik- 
ing ways in which a quantum system can show deviation 
from ergodicity at the single-channel level even though 
the classical dynamics is completely ergodic. Wavefunc- 
tion intensities near a short unstable periodic orbit follow 
a distribution far from that predicted by random matrix 
theory (RMT), with some wavefunctions having much 
more intensity and other much less than would be pre- 
dicted based on gaussian random fluctuations. The phe- 
nomenon is at first glance paradoxical, because the long- 
time (and indeed stationary) quantum behavior retains 
a memory of the short-time classical motion, a memory 
that is completely absent in the long time classical dy- 
namics of a chaotic system. Scarring has been observed 
experimentally in a wide variety of systems, including mi- 
crowave cavities semiconductor structures M, and 
the hydrogen atom in a magnetic field [Bp). 
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A theory of scarring based on the linearized evolution 
of gaussian wavepackets was firstprovided in later 
theoretical work by Bogomolny J7J in coordinate space 
and Berry || in Wigner phase space followed. These 
made predictions about the average intensity on a clas- 
sical periodic orbit of states in a given energy band; 
however, because of the energy smoothing involved, no 
predictions were possible about the statistical properties 
of individual peak heights in the local density of states. 
More recently a nonlinear theory was developed Q which 
made it possible to predict the statistical properties of in- 
dividual wavefunctions, in the semiclassical limit. A ho- 
moclinic orbit analysis showed that long-time return am- 
plitude to the vicinity of a periodic orbit bore the imprint 
of the short-time linearized classical dynamics around the 
periodic orbit. This leads to a natural separation of scar- 
ring intensity into a classical short-time component and 
a random long-time component, as suggested already in 

In predictions were made about the distribution 
of wavefunction intensities on a periodic orbit and at a 
generic point in phase space. The full distribution of in- 
tensities, which includes samples taken over all of phase 
space, has a long tail (compared to the Porter-Thomas 
prediction of RMT), dominated by the effect of the least 
unstable periodic orbit. The functional form of this tail is 
given in the semiclassical (high-energy) limit very simply 
in terms of the stability exponent of this least unstable 
orbit, as long as an optimally-oriented test basis is cho- 
sen. Furthermore, upon ensemble averaging a power-law 
intensity distribution tail is obtained, in sharp contrast 
with the exponential tail predicted by RMT. This result 
is also to be contrasted with the log-normal intensity dis- 
tribution tail which obtains in diffusive systems p2|^3[ . 
Thus, although RMT is accepted as the zeroth-order ap- 
proximation for both chaotic and disordered quantum 
systems, deviations from RMT predictions can be quali- 
tatively different in the two cases, providing an impetus 
for the present research. 

The numerically tested quantitative predictions in 
concerned the local densities of states of closed 
systems. Experimentally, certain properties of open sys- 
tems, such as resonance widths and conductances may be 
more amenable to experimental verification JTif . Much 
work has been done here within the context of RMT (jl5| ; 
deviations from RMT in disordered systems have also 
been considered In the present work, we address 

the distribution of decay lifetimes in a leaky chaotic sys- 
tem, and the probability to remain in such a system [17] ] 
as a function of time and the location of the "leak" (open 
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channel) . The distribution in phase space of the remain- 
ing probability density at long times is considered, as well 
as the dependence of the probability to remain on where 
the particle was first "injected" into the system. These 
last two questions bring us into contact with wavefunc- 
tion correlations and transport in chaotic systems, which 
(as we show) can be very different from RMT expecta- 
tions where periodic orbits are involved. All the quanti- 
tative predictions which follow are tested numerically. A 
study of conductance properties in two-lead chaotic sys- 
tems, including mean conductance, conductance fluctua- 
tions, distribution of peak heights, and peak correlations 
(and how all these depend on the placement of one or 
both leads in relation to the classical orbits) is forthcom- 
ing jl|. 

II. CLASSICAL AND QUANTUM WEAKLY 
OPEN CHAOTIC SYSTEMS 

A. Classical behavior 

We begin by considering a small opening in a classi- 
cally chaotic system, which allows a particle to escape 
from the system. We will often use language suggest- 
ing that the "opening" is defined in position space, as it 
would often be, for example, in a mesoscopic experiment. 
However, the formalism considered here is much more 
general: all that is required is that the opening be local- 
ized in the classical phase space; escape routes that are 
defined exclusively in terms of position or momentum are 
special cases of this. A simple example of a momentum 
space opening is a potential barrier that allows particles 
to leave only if their momentum is directed almost nor- 
mal to the wall plfl . An opening having the shape of a 
phase-space gaussian naturally occurs when one consid- 
ers tunneling out of a metastable chaotic well formed by 
a continuous potential fl9|| . Now we can imagine forming 
a mesh in classical phase space with each cell the size 
of the opening; because the classical dynamics (in the 
closed system) is chaotic, probability density starting in 
one such cell will soon be evenly distributed over all the 
available cells. The time for this to happen is logarithmic 
in the size of the opening w. 

T mix - i| logw| . (1) 
A 

Here A is the Lyapunov exponent of the classical dynam- 
ics (the mean rate of chaotic divergence of classical or- 
bits), and the total size of phase space, in terms of which 
w must be measured, has been set to unity. On the other 
hand, the escape time from the system is inversely pro- 
portional to the leak size w, so a small value of w will 
cause complete mixing of the remaining probability to 
take place on a time scale much shorter than the scale 
on which probability is leaking out. One obvious conse- 
quence is that the probability to remain in the open clas- 
sical system follows an exponential law. This behavior is, 



of course, independent of the position of the leak. Also, 
density is constantly redistributing itself, so that the re- 
maining probability density remains evenly distributed 
on the scale of our mesh, except in a corridor of length 
scaling as 1/ A leading forward in time from the position 
of the opening. Notice that the width of this corridor 
where the probability to remain is suppressed scales as 
the size of the opening. Thus, this corridor has no ef- 
fect on the quantum behavior when the opening is small 
compared to fi. Finally, even if the initial probability is 
not evenly distributed over the entire phase space, the 
long-time behavior is unaffected (as long as the bulk of 
the probability is not initially placed in a corridor simi- 
lar to the one described above, but leading backwards in 
time from the opening) . 

In contrast to these results, we will find in the quantum 
case that the probability to remain in the system at long 
times depends strongly on whether the opening is located 
on a classical (unstable) periodic orbit, even though the 
initial probability density is evenly distributed. Again, 
we see that long-time quantum behavior retains a bet- 
ter memory of short-time classical dynamics than does 
the long-time classical behavior. Also, we will see that 
given a leak placed on a periodic orbit, the remaining 
probability distribution at long times can be strongly af- 
fected not just on the periodic orbit itself, but also on the 
other short periodic orbits of the system. Enhancement 
or suppression can be observed depending on the energy 
range considered and on the classical actions of the or- 
bits in question. Similarly, the probability to remain at 
long times will be affected if the original probability is 
injected on an unstable periodic orbit different from the 
one where the opening is located. All this is true even 
though the decay is taking place on a time scale much 
longer than any other time scale in the problem (the pe- 
riod of the short orbit, the mixing time, and also the 
Heisenberg time, i.e., the inverse of the mean level spac- 
ing). 

B. Quantum mechanics and RMT 

Let the quantization of our classical system be given 
by an A— dimensional Hilbert space (N is the number 
of Planck-sized cells in the classical phase space), with 
dynamics in the closed system given by the Hamiltonian 
Hq. If the opening is very small (less than one open 
channel, so that the resonances are non-overlapping), we 
can write an effective Hamiltonian for the open system: 

H = H -i^\a){a\, (2) 

where \a) is a quantum channel associated with the open- 
ing, and r is the decay rate in that channel (taken to be 
small). | a) could be a gaussian wavepacket enclosing the 
hole, or a position or momentum state. It is important 
to note here that the opening is small compared to h 
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(less than half a wavelength if in position space). One 
can of course consider large openings, or ones which are 
not thus localized to a single channels; these possibili- 
ties are considered towards the end of the present paper. 
We emphasize, however, that the phenomenon discussed 
here is already present in its full form for the tiniest sin- 
gle opening, without the complications that arise in the 
more general case. 

For a small hole, the main effect of the opening on 
a wavefunction |\& n ) of the closed system is that it ac- 
quires a decay width proportional to the intensity of the 
wavefunction at the opening: 



r, 



r|(*n|a)| 5 



(3) 



If the intensities x n = N\(^ n \a)\ 2 follow a chi-squared 
distribution, as in RMT, we have probability P(x) = 
-j= — exp (— x/2) for real overlaps (^ n \a), and P(x) — 
cxp (— x) for complex overlaps. Consider the complex 
case. Because mixing between the states of the closed 
system can be neglected in the small V regime, the total 
probability to remain in the system is given by a sum 
over these states: 



jV-1 



dxP{x)e 



-xTt/N 



dxe- x e- xTt ' N 



1 



i + vt/N 



(4) 



[Remember that N is the total number of states in the 
system; the classical decay rate is given by T c i = T/N 
because only one channel has the possibility to decay] 
We see that at short times (t <§; T^ 1 ), the probability to 
remain in the system is P iem (t) ~ 1 — T c \t, as expected, 
while at long times we have the asymptotic behavior 



Pr, 



(5) 



In the case of M independent weakly open channels, i.e. 



M-l p (i ) 

T 



i=0 

the classical decay rate is given by 

M-l 

N 



M-l 

r i - — V r w 



i=0 



and the RMT probability to remain is 

M-l 

*W = n TTr^ 



i=0 



*)t/N 



(6) 



(7) 



(8) 



Taking M — > oo while keeping the total decay rate T c i 
constant, exponential decay consistent with the classical 
prediction is obtained. On the other hand, fixing the 
number of channels M and taking t — > oo, we observe 
the power-law behavior 



.(*) 



(N/t) 



M 



(9) 



The case of real overlaps (W„|a) follows similarly: each 
real random overlap counts as half of a complex one, 
so there P rem {t) <~ t~ M l 2 . In the literature one often 
considers the distribution of delay times for scattering 
off of the system in question: there one must include 
the probability of populating a given resonance in the 
first place, which of course is proportional to r„. This 
leads to an extra factor of t in the denominator, giving 



P. 



delay 



(t) 



t 



-M- 



1 for complex overlaps and fdeiay(i) 



£-M/2-i £ or rea j over i a p S> j n our caS e ; we imagine the 
system to be populated first, before the lead is opened 
up, and thus no extra power of t is present. 



III. EFFECT OF PERIODIC ORBITS 

A. Probability to remain 

We now go beyond RMT to consider the effect of real 
dynamics on the quantum probability to remain in a clas- 
sically chaotic system. Take the channel \a) to be on 
or near an (unstable) periodic orbit of instability expo- 
nent A. The smoothed local density of states at \a) is 
obtained by fourier transforming its short-time autocor- 
relation function, which is easily obtained by linearizing 
the classical equations of motion near the unstable or- 
bit ^,^,0. Thus, for example, if the periodic orbit in 
question is a fixed point of a discrete time map, and \a) 
is a gaussian wavepacket optimally aligned along the sta- 
ble and unstable manifolds of the orbit, the short-time 
autocorrelation function is given by 



A U n(t) = (a\a(t)) = 



-i(pt 



V cosh Xt 



(10) 



Here —(f) is a phase associated with one iteration of the 
orbit: it is given by the classical action in units of h, 
plus a Maslov phase as appropriate. The subscript "lin" 
indicates that the expression is obtained within the lin- 
earized classical approximation; it is valid on time scales 
short compared to the mixing time T m ; x ~ =| log ?i | . 

A more general form of Eq. [lO] applies for a non- 
optimally oriented wavepacket (e.g. a position or mo- 
mentum channel could be non-optimal depending on the 
direction of the invariant manifolds at the periodic point), 
and also for a channel not exactly centered on a periodic 
point |l^,|2(|. In particular, for a wavepacket centered on 
the periodic orbit but not optimally oriented with respect 
to its invariant manifolds, the form above becomes 
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Au a (t) = (a\a(t)} 



-i(pt 



\/cosh Xt + iQ sinh Xt 



(ii) 



In Eq. [ll], Q is a non-optimality parameter: in a coor- 
dinate system where the stable and unstable manifolds 
are orthogonal, Q is a function of the angle between the 
orientation of the phase-space gaussian (at some fixed 
eccentricity) and either of these two directions. Alter- 
natively, if the wavepacket \a) is fixed to have a circular 
shape in phase space (i.e. to have equal and uncorre- 
cted uncertainties in q and p), Q becomes a function 
of the non-orthogonality between the stable and unsta- 
ble manifolds. In any case, as long as Q is not very 
large, the qualitative behavior is not much changed, and 
although analytic results are less easy to obtain for non- 
zero Q, quantitative predictions can be readily produced 
for comparison with any experimental or numerical data. 

The key point for our purposes here is that for a small 
exponent A, the autocorrelation function remains large 
for the first 0(A _1 ) iterations of the orbit, and the local 
density of states has a short-time envelope 



S lin (E) = ^2e iEt A nn (t) 



(12) 



of width scaling as A and height scaling as A -1 [also see 
Fig. [l] below] . 

Nonlinear recurrences on time scales beyond the mix- 
ing time (associated with orbits homoclinic to the orig- 
inal periodic orbit) lead to fluctuations multiplying this 
spectral envelope, eventually producing a line spectrum 



S(E)=J2\{*n\a)\ 2 5(E-E n ). 



(13) 



The line heights x n = N\(^> n \a}\ 2 are distributed in each 
energy region as a chi-squared distribution with mean 
SUE) Q: 



P(x) 



1 



,-x/S Un (E) 



(14) 



Thus, the distribution of decay widths can be strongly 
energy-dependent; in particular, the probability to re- 
main in the system at long times is now given by 

/>oo 

P rem (i) = / dxP(x)e- xTt ' N 
Jo 

1 



1 + S lin (E)Tt/N 

1 1 



Slin{E) T c it 



(15) 



if initially only states with energy around E are popu- 
lated. 

The scarred states (those with energy close to satis- 
fying the EBK quantization condition) have Sy ln (E) > 
1 and thus decay much faster than the antiscarred 
states [EOl, which are far from satisfying EBK and thus 



have S\i n (E) < 1. Let us examine more closely these two 
distinct energy regimes. Near the quantization energy 
E = (p, the smoothed density of states Sn n (E) has its 
peak; its height scales inversely with A for small A (|] : 



Slin(E — 



c/A, 



(16) 



where c = 5.24 is a numerical constant [fLl| . The width [] 
of this peak in S\i n (E) scales linearly with A for small A, 
and all of the anomalously enhanced wavefunction inten- 
sities come from this energy region, as was observed and 
confirmed numerically in jll|. In the open system, these 
states produce an excess of large resonance decay widths 
and decay faster (by a factor of C^A" 1 )) than would be 
predicted by RMT. 

Because our focus here is on the long-time behavior of 
weakly open systems, we are more interested in the (com- 
plementary) suppression of the smoothed local density of 
states far from the resonance energy. Again, we consider 
the strong scarring regime, where A <§; 1: then the linear 
spectrum falls off exponentially far away from the peak, 



Slin (E) » ^ e --|s-0l/2A 



(17) 



for \E— 4>\ ^> A. Within O(A) of the optimal anti-scarring 
energy, E = <p + tt, the spectrum deviates from the expo- 
nential law and smoothly approaches the value 



S lia (E = <f> + n) 



47T 



-tt 2 /2A 



(18) 



at the minimum. The region within O(A) of E = <p + tt is 
thus responsible for producing the smallest wavefunction 
intensities, and the narrowest resonances in the corre- 
sponding open system. This excess of exponentially small 
decay rates is as dramatic a signature of the underlying 
classical behavior as the long wavefunction intensity tails 
found in pl| . As we will observe in the next section, the 
antiscarring effect on the long-time behavior of open sys- 
tems can be very striking even for moderate exponents A 
(e.g. A w 1), as long as the lead is optimally placed with 
respect to the periodic orbit. 

Smoothed local densities of states Su n (E) on periodic 
orbits of instability exponents A = 0.20 (solid curve) and 
A = 0.15 (dashed curve) are shown graphically in Fig. [I]. 
The figure can be viewed as representing the average 



We note that our presentation here is in the context of a 
discrete-time map; thus E is a dimensionless quasienergy that 
takes values in the interval [0, 2n]. For a real continuous-time 
system with a periodic orbit of period Tp, it is of course the 
quantity ETp/h that must be compared with the dimension- 
less number A. Also, the smoothed density of states will then 
have an infinite sequence of peaks, each centered on an energy 
satisfying the EBK quantization condition [[|. The ratio of 
each peak width to the spacing between peaks scales as A. 
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wavefunction intensity in a closed system as a function 
of energy, or the mean resonance width (in units of T c i) 
at that energy in the weakly open system. The phase 
4> = has been chosen so as to make E = the EBK 
energy at which maximum scarring occurs (Eq. 16). A 



half-log scale is used to emphasize the exponential falloff 
in average resonance width between E = and E = tt 
(Eq. [I?]), and the minimum near the anti-EBK energy 
E = tt (Eq. |l^). For reference, the smoothed local den- 
sity of states in RMT (applicable when the lead is not in 
the vicinity of any short periodic orbit) is displayed as a 
dotted line in the figure. 

We now consider the energy-averaged probability to re- 
main in the open system: this will be the quantity stud- 
ied in detail numerically in the next section, where the 
model system is a (non-energy conserving) discrete-time 
kicked map. [For an energy-conserving system, varying 
the strength of a weak magnetic field and thus sweeping 
through different values of the phase cj) would produce 
the same result.] Again, because the perturbation in- 
duced by opening up the system is small, there is little 
mixing among states of different energy. Thus the total 
probability to remain is obtained simply by averaging the 
probabilities at the different energies. From Eq. [l5| we see 
that at short times, the classical behavior is recovered: 



Prem = 1- < ftin > Tt/N = 1 - T cl t . 



(19) 



as < Sl in >= Ai in (0) 



1 by normalization. 



{a\a) 

Thus at short times, t <C r cl 1 , the faster-decaying scarred 
states and slower-decaying antiscarred states always can- 
cel exactly and no sign of the quantum signatures of the 
underlying classical dynamics can be observed. On the 
other hand, at long times, i.e. for t ^> [Siin(-Emin)r c i] _1 , 
we obtain the very different behavior 



< S,7 



> 



(20) 



Here E m i n is the energy at which the smoothed spectrum 
has its minimum; for an optimally placed lead \a) this 
energy is exactly tt out of phase with the EBK energy 
(j>, as discussed above (Eq. ^8|). < > is the inverse 
of the smoothed density of states at \a), averaged over 
energy (or weak magnetic field). 

As < S]ja >= 1 by definition, any fluctuations in the 
smoothed spectrum resulting from short-time recurrences 
will cause < > to be greater than one, resulting 

in an enhanced probability to remain at long times. In 
particular, for an optimally placed lead (corresponding 
to Eq. |l^), let us consider the strong scarring regime of 
small A. This gives the exponentially large enhancement 



7r 2 /2A 



(21) 



This long-time behavior is completely dominated by the 
most antiscarred states, i.e. those with energy within 
0(A) of E min = <f> + 7r (Eq. 111). For a non-optimally 



placed lead, with not too large non-optimality parameter 
Q (see Eq. [ll] and discussion following) we find empiri- 
cally a similar exponential enhancement of the typical 
decay time: 



(ir 2 /2-bQ)/\ 



(22) 



where b = 1.1 is a numerical constant. 

If the state \a) defining the phase-space location of the 
opening is centered off of the periodic orbit, but within Ti 
of the orbit, one still has fluctuations in the linear density 
of states and consequently an enhancement in the prob- 
ability to remain at long times. An analytic form for 
the linear autocorrelation function in such a case can be 
found in p0| . For a circular minimum-uncertainty phase- 
space opening centered a distance 5 away from a periodic 
orbit with small exponent A, the energy-averaged value 
< S^} > scales as 



< ^lin 



)2 e (ir 2 /2-dS/Vh)/\ 



(23) 



where d is yet another numerical constant. S can be a 
displacement along either the stable or unstable direc- 
tion away from the orbit. Thus, deviations from RMT 
behavior are observed in an area scaling as Ti surrounding 
the periodic orbit. Maximum enhancement of (i.e. 

enhancement of order A e" ' ) occurs for S < O(XVh), 
corresponding to a phase-space area scaling as A 2 h sur- 
rounding the orbit. Thus, if we consider the long-time 
probability to remain in the system averaged over all pos- 
sible positions of the lead, we obtain 



p — 

J rem — 



i + o(fiAy /2A ) 



(24) 



[The correction to RMT is obtained by multiplying the 
maximum obtainable enhancement by the size of the 
phase-space region where such enhancement occurs.] In 
principle, contributions from all the periodic orbits need 
to be added, however, if orbits with small A exist, they 
will clearly dominate any such sum. The result is that 
at finite energy, exponentially large (in 1/A) deviations 
from RMT are found even in the phase-space averaged 
analysis. In the h — > limit of any given classical system, 
the RMT behavior is recovered because the chance of a 
lead being fo und o n the short periodic orbit goes to zero. 

In Section IV B , we present theoretical predictions and 
numerical data measuring the probability to remain in 
the system at long times as a function of the location of 
the opening. 



B. Probability density at long times and dependence 
on initial conditions 

Up until now we have been focusing on the distribu- 
tion of resonance widths and on the total probability to 
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remain in the system starting from a uniform initial state, 
all as a function of the location of the lead \a). In other 
words, while changing the location of the opening, we 
have always been tracing over the initial and final states 
of the system. We now proceed to address two related 
questions, both of which require us to consider transport 
properties and wavefunction intensity correlations within 
the system. 

First, still taking the initial filling to be uniform, and 
fixing the location of the opening to be on a periodic 
orbit, it is natural to ask how the long-time probabil- 
ity density remaining in the system distributes itself over 
phase space. Classically, we expect the remaining prob- 
ability always to be redistributing itself on a time scale 
short compared to the decay time, and thus to be uni- 
form except in a very narrow corridor encompassing the 
unstable orbit. The width of the corridor scales as the 
size of the opening. In RMT, of course, the remaining 
probability is also completely uniform except at phase 
space locations having non-zero overlap with \a). In con- 
trast, we find that in the real quantum system, the re- 
maining probability density is strongly suppressed in a 
corridor of size h around the orbit, much wider than the 
size of the lead. Even more interesting is the fact that the 
probability density can be either relatively enhanced or 
suppressed along the other unstable orbits of the system, 
depending on the classical actions associated with these 
orbits. 

Before proceeding, we mention a closely related prob- 
lem, which can be thought of as a time-reversed version 
of the one stated above. Instead of initially filling the sys- 
tem with a uniform density, we inject probability in some 
known initial state and look at the probability to remain 
after a long time as a function of this initial state. This 
state, which we call \b), should be classically well-defined, 
i.e. it can be a phase-space gaussian, or a position or mo- 
mentum state, as discussed above. 

The two problems are in general distinct: if H is the 
non-Hermitian quantum Hamiltonian, the first involves 
the quantity (b\e~ lH t e lHt \b) 1 while the second measures 
(b\e lH t e~ lHt \b). However, when T is very small (in the 
regime of non-overlapping resonances), H is nearly nor- 
mal, the distinction between left and right eigenstates 
vanishes, and the two quantities both converge to the 
eigenstate sum 



E 



l*n> 



-r„t 



(25) 



For | b) not on the periodic orbit containing the lead \a), 
the quantity |(&|'I' n )| 2 is independent of r„ ~ |(a|^„)| 2 , 
and follows its own chi-squared distribution with mean 
scaling as S# n (E). Here S^ in (E) is the fourier transform 
of the linearized (short-time) autocorrelation function of 
the test state \b); it is to be distinguished from Su n (E) = 
Sy in (E), the smoothed local density of states at the lead. 
We easily obtain, at energy E, 



l + S% n (E)rt/N 
(E) 1 



°1i 



SL(E) r cl t 



(26) 



Averaging over E, we obtain the ratio of the remain- 
ing probability density at \b) to the average remaining 
density at long times: 



/SL > 



(27) 



We see that this ratio goes to unity if S^ in has no energy 
dependence, i.e. if \b) does not lie on a short periodic or- 
bit. If the position \a) of the lead itself does not lie on a 
periodic orbit, the remaining density profile will of course 
be flat over all states \b). However, if both \a) and \b) lie 
on periodic orbits, the probability to be found at \b) can 
be either suppressed or enhanced, depending on whether 
the energy envelopes Sy ln and S^ in are in or out of phase in 
the energy range being averaged over. For simplicity, let 
us consider an example where the periods, and instabil- 
ity exponents of the two orbits are equal. Then the two 
smoothed energy envelopes are identical, up to a relative 
phase shift (the difference between <p a and <pb), which can 
be adjusted by varying a magnetic flux enclosed by one 
of the orbits. If the two are exactly in phase, S 1 ^ = St , 
then the ratio in Eq. |27| reduces to 1/ < >, which, 
we recall from our previous discussion, is a quantity ex- 
ponentially small in the instability exponent A. Thus, 
the remaining probability very strongly avoids the orbit 
on which \b) is located. Another way of expressing this 
result is that the total probability to remain in the sys- 
tem at long times is exponentially suppressed if the initial 
state is located on an orbit which is "in phase" with the 
orbit on which the opening is located. 

The suppression of probability density given by Eq. |2?] 
is of course a pure quantum interference phenomenon; 
it has no analogue in the classical dynamics of open 
systems. It is also fundamentally a long-time effect as 
there is in general no short path leading from \a) to \b) 
which could give rise to such intensity correlations. How- 
ever, despite being intrinsically long-time and quantum, 
the phenomenon can be understood only in terms of the 
short-time, classical dynamics near each of the two un- 
stable periodic orbits. This demonstrates once again the 
power of semiclassical techniques for understanding long- 
time quantum behavior. 

In the opposite extreme case, where the two orbits are 
out of phase_exactly by tt [S? in (E) = S? m (E + tt)], the 



ratio in Eq. 27 is dominated by the region of the enve- 



lope where Sf~ is maximized and Su n minimized. [This 
is an energy region in which the wavefunctions tend to be 
scarred near \b) and antiscarred near |o}.] The relative 
intensity enhancement at |6) then scales with the height 
of the peak in Sf in , i.e. as A -1 ^> 1. So a large enhance- 
ment of the remaining probability is found on orbits out 
of phase with the one on which the opening is located. 
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We need to consider also the case where states \a) and 
1 6) are found on the same orbit [the same reasoning ap- 
plies if | a) and \b) are on distinct orbits that are related by 
a symmetry transformation]. This corresponds to mea- 
suring the remaining probability along the orbit on which 
the lead is located (or alternatively to launching the ini- 
tial probability along this orbit). First, consider the case 
where \a) and \b) are exactly related by time evolution 
in the closed system. Then the two local densities of 
states are identical, i.e. |(a|^„)| 2 = | (&| ^I/n) | 2 for each 
ri. It is easy to see from Eq. p5| that P Iem in this case 
decays at long times as 1/t 2 instead of the usual 1/t be- 
havior. This is easy to understand intuitively: the very 
long-lived resonances which survive at long times have 
very little amplitude at \b). More generally, let us con- 
sider | a) and |6) lying on the same orbit but not exact 
time-iterates of one another. This is possible even if \a) 
and | b) are both optimal (in the sense of having Q = 0, 
see Eq. [ll]). Thus, the iterates of \a) may have width 
<7oe A ™ along the unstable manifold as they pass through 
that point on the orbit on which \b) is centered p2fl . If 
we choose a width for \b) which does not correspond to 
any integer n, then \b) is not any exact time-iterate of 
| a). However, for some time i we may still write 



\b)=a\a{t))+ 1 \c), 



(28) 



where \a\ 2 + \j\ 2 = 1. Then the local density of states 
at \b) separates naturally into two parts: one of weight 
\a\ 2 which is exactly equal to the density of states at 
the opening \a), and another of weight 1 — \a\ 2 which is 
statistically independent of the former but has the same 
linear energy envelope. The first, as we just saw, gives 
a contribution to P r b cm which scales as 1 /t 2 and thus can 
be ignored at long times. The second behaves just as 
if 1 6) were located on a different orbit having the same 
linear envelope. Thus for \a) and \b) on the same orbit we 
obtain the same exponential suppression factor as before 
(Eq. |27| ), times the extra suppression factor 1— \ a\ 2 . This 
latter factor also becomes very small for small A [ p2[ , 
as any wavepacket optimally placed on a periodic orbit 
comes ever closer to being an exact time-iterate of any 
other such wavepacket on the same orbit. 

We note again that this effect is purely quantum- 
mechanical, based though it is on short-time semiclassical 
analysis. Classically, only a tiny fraction of the proba- 
bility distribution corresponding to \b) would leak out 
through the hole at \a) before the density escapes from 
the periodic orbit and proceeds to distribute itself evenly 
over the entire accessible phase space. 



IV. NUMERICAL TESTS 

A. The model 

We now proceed to test numerically the various results 
obtained analytically in the previous section. What is 



required is a large ensemble of chaotic systems with each 
realization having a short unstable periodic orbit of the 
same instability exponent. For this purpose we consider 
kicked maps on the toroidal phase space [—1/2,1/2] x 
[—1/2, 1/2]. The classical dynamics for one time step is 
given by 



p — > p = p + mq — V (q) mod 1 
q —y q = q + np + T (p) mod 1 . 



(29) 



This dynamics can be obtained from the stroboscopic 
discretization of a kicked system 23 1 with a kick potential 
-|mg 2 + V(q) applied once every time step and a free 
evolution governed by the kinetic term \np 2 + T(p). m 
and n are arbitrary integers, while V, T are periodic 
functions of position and momentum, respectively. The 
system can also be thought of as a perturbation of the 
linear system (cat map) p4j 



p — ► p = p + mq mod 1 

q — ► q = np + (mn + l)q mod 1 . 



(30) 



For given positive integers m, n, we choose the func- 
tions V and T such that m — V"(q) > for all q and 
similarly n + T"(p) > for all p. Then the system is 
strictly chaotic and looks everywhere locally like an in- 
verted harmonic oscillator. 

The quantization of such systems is well-studied in the 
literature p3| . h must be chosen such that N — Xjl-nti, 
the number of h— sized cells in the classical phase space, is 
an integer (N must be even to preserve the periodicity of 
the quadratic terms in the potential and kinetic energy) . 
For doubly-periodic boundary conditions, the quantum 
N— dimensional Hubert space is spanned by the position 
basis where qi — i/N and i = 0...N — 1. The 
momentum-space basis is given similarly by \pj), where 
Pj = j/N and j = . . . N — 1; and the two bases are 
related by a discrete fourier transform. The quantum 
dynamics is then given by a unitary N x N matrix 



U — exp 
• exp 



- / ( - 2 „ir ■ l ip ) 
/ ( -mq 2 - V(q) } ■ h 



(31) 



Each factor is evaluated in the appropriate basis, and 
an implicit forward and backward fourier transform has 
been performed. 

We may now perturb this unitary dynamics to allow 
for a small decay rate in channel a: 



W=(l-^\a)(a\\r. 



(32) 



Eq. |32j is of course the discrete-time version of the 
continuous-time dynamics given in Eq. ^| above. Since 
we are working in the regime r C 1, where the decay 
rate per time step is small, the discretization of the de- 
cay process will not affect the long-time behavior of the 
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system. The decay channel \a) can in principle represent 
any vector in the Hilbert space; however, for the decay 
to correspond to a classical escape route, \a) should be 
a phase-space localized state, such as a position or mo- 
mentum state. We will find it convenient to let \a) be a 
circular phase-space gaussian in the coordinates q, p. 

We now need to construct an ensemble of such systems, 
all having the same behavior in the vicinity of a short 
unstable periodic orbit. For this purpose, we set m = 
n = 1, and let the potential V and kinetic term T be odd 
functions of their respective arguments: 



V(q) = ^2[K r sin(27rrg) - 2tt K r rq] 

r=l 
3 

T{p) = sin(27rrp) - 2n;K' r rp] 



(33) 
(34) 



r=l 



It is easy to see that the equations of motion (Eq. |29| ) 
then have a fixed point at the origin with Jacobian matrix 

' ' [the system can be thought of as a perturba- 



J = 



1 2 
tion of the 



cat map, with perturbation vanishing 



1 1 
L 1 2 

near iq,p) = (0,0)]. The instability exponent is given by 
A = cosfT 1 (^Trj) = 0.96. We may now choose each of 
the coefficients K r , K' r from a uniform distribution over 

0.3 0.3 



the interval 



One easily sees that each 



' (2-rrr) 2 ' (27rr-) 2 

system in this ensemble satisfies everywhere the condi- 
tion 1 — V" > 0, 1 + T" > mentioned above, which is 
sufficient to ensure hard chaos. 



B. Probability to remain 

The ensemble-averaged probability to remain in the 
system after t time steps is now computed for various 
positions of the exit channel \a). For simplicity, we choose 
| a) to be a circular phase-space gaussian 



a(q) 



-(Q-Qo) /2h+ip (q-q Q )/h 



(35) 



centered on (qo,Po) and having width yh in both the 
q and p directions. Because the Jacobian J at the fixed 
point (0, 0) is symmetric (and the stable and unstable di- 
rections are thus locally orthogonal), such a wavepacket, 
when centered at (0,0), is optimal in the sense of hav- 
ing the slowly decaying short-time autocorrelation func- 
tion of Eq. [HJ [A more general gaussian wavepacket 
(including position or momentum states as extreme lim- 
its) centered on the periodic orbit may have a non-zero 
parameter Q (see Eq. [ll]), leading to a less sharp linear 
spectral envelope and less strong scarring and antiscar- 
ring effects. The qualitative behavior would, however, 
remain unchanged.] 

In Fig. |^, the probability to remain in the system as 
a function of the scaled time t' = T c \t is first plotted 



(using plusses) for a generic leak location. The data was 
collected for systems of size N = 120 and decay param- 
eter T = 0.1. The results agree well with the the RMT 
prediction P rom (0 = 1/(1 + (dashed curve). For com- 
parison, the classical probability to remain, exp(— t'), is 
plotted as a dotted curve. Next, we place the opening on 
the periodic orbit at the origin of phase space, and obtain 
the rather different behavior, with an enhanced long-time 
tail (squares). The asymptotic form is well reproduced 
by the scar theory prediction, -P rcm (0 ~< ^lln > /O 
which is shown in Fig. |^ as a solid line. For the instabil- 
ity exponent A = 0.96, we observe a long-time probability 
enhancement factor < S^ 1 >= 11.04. Of course, bigger 
enhancement factors can be observed for less unstable 
orbits, as we will see below. 

First, we examine more carefully the probability to 
remain at long times as a function of the position of 
the lead. In Fig. ||a is plotted the total probability to 
remain in the system at time t' — 2 • 10 3 , for vari- 
ous locations of the opening (all for N = 30). These 
possible locations are located on a 40 x 40 grid filling 
the middle l/9th section of the total phase space (i.e. 
iq,p) G [-1/6,1/6] x [-1/6,1/6]). The bright spot at 
the center of the figure represents the enhanced prob- 
ability to remain if the opening is located exactly on 
the periodic orbit. As we see from the figure, the an- 
tiscarring effect falls off quite quickly as the opening is 
moved away from the periodic orbit (in fact, the size of 
the bright spot scales as h; see Eq. |2^ and discussion 
following). In Fig. ||b is plotted the theoretical quantity 
< S^ n >, as computed using the linearized equations of 
motion (Eq. |3(]) around the periodic orbit. This is ob- 
served to be in good agreement with the data. Of course 
the linearized equations of motion only hold near the pe- 
riodic orbit itself, and do not correctly describe classical 
motion in other regions of phase space. However, in our 
case the short periodic orbit at the origin clearly domi- 
nates the data. If the classical system contained several 
not very unstable orbits (see next subsection for an exam- 
ple), several bright spots would appear in the plot, and 
each could then be well reproduced using the linearized 
classical dynamics around the appropriate orbit. 

An important feature to notice in Fig. ||is that P rom (t) 
at long times depends not only on the distance of the lead 
from the periodic orbit but also on the direction. Greater 
enhancement is observed if the lead is placed along cither 
the stable or the unstable manifold of the orbit (the two 
"diagonals" of the would-be "square"). 

We now consider how the probability to remain at long 
times depends on the instability exponent A of the orbit 
on which the lead is located. As we discussed in the pre- 
vious section, as A gets small, resonances exponentially 
narrow in A should appear at the antiscarring energies, 
and the total probability to remain at long times is en- 
hanced by a factor exponentially large in A (Eq. |l]). To 
see this effect, we modify the classical dynamics of Eq. ^ 
by adding an additional term to the potential and kinetic 
functions: 
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V(q) = ...- 



A 



cos(27r<7) 



(2^) 2 

T{p) = ■■■ + C0S(27T P ) . 



(36) 



The same value of A should be used in the potential and 
kinetic terms to preserve the symmetry of the Jacobian. 
The Jacobian matrix of the dynamics near the periodic 
orbit at (0, 0) is then given by 



1 1- A 

1-A 1 + (1-A) 2 



(37) 



For positive A, the trace of the Jacobian decreases and 
the orbit becomes less unstable. Using A = 0.0, 0.1, 0.2, 
0.3, and 0.4, we obtain exponents 0.96, 0.87, 0.78, 0.68, 
and 0.59, respectively. 

In Fig. H the long-time enhancement factor of P rom (t) 
over its RMT value is plotted as a function of the ex- 
ponent A, using plusses for N — 120 and squares for 
N = 240. The theoretical prediction < > is shown 
as a solid curve. The data consistently falls below the 
theoretical prediction, with the disagreement becoming 
more pronounced at the smaller values of A. The reason 
is primarily a finite-size effect: the analytical calcula- 
tions are all carried out under the assumption that the 
mean level spacing is much smaller than any scale over 
which the linear energy envelope changes significantly. 
Then the linear envelope is roughly constant on the scale 
at which individual resonances emerge, and their behav- 
ior can be treated statistically. Thus, the discrepancy 
becomes more noticeable as A — * for fixed AT, as the 
structures in Su n become more comparable to the mean 
level spacing. Indeed, we see that the N = 240 data is 
consistently closer than the N = 120 data to the theory, 
which strictly applies only in the semiclassical N — > oo 
limit. 

We observe the exponential increase in the enhance- 
ment factor as A decreases; indeed the very moderate ex- 
ponent A = 0.59 produces an enhancement factor of well 
over 100 in the long-time probability to remain. A = 0.1 
would in theory produce an average long-time enhance- 
ment of 1.2 • 10 13 , provided we were able to go to a large 
enough system, wait for a long enough time, and collect 
enough statistics to observe it [from Eq. [l^ we see that 
the asymptotic 1/Sn n t' form holds only for t' ^> S^ n ]. 



C. Phase space distribution of long-time probability 
density 

We now consider predictions concerning the enhance- 
ment or suppression of the long-time probability near or- 
bits other than the one on which the lead is located [see 



discussion in Section IIIB]. To test these predictions we 
need an ensemble of systems all having two periodic or- 
bits in common, and the ability to vary the action phase 
difference between them. In our example the two orbits 



the Jacobian J 



have the same local dynamics in their respective neigh- 
borhoods, though this of course is not necessary to pro- 
duce the desired effect. 

We work on the phase space (q,p) € [—1/4,3/4] x 
[—1/2, 1/2], set m = n = 2 in the equations of motion 
(Eq. |2^), and impose the constraint K\ = — 3K^ on the 
kick potential (see Eq. |33| ). This condition ensures the 
presence of a fixed point at (1/2,0) in addition to the 
usual one at (0,0) on which we have been focusing so far. 
The linearized dynamics around each orbit is given by 

1 2 1 

2 g , and the exponent per period 

is A = 1.76. Other orbits are of course present, but they 
change with the coefficients K r , K' r , and so their effects 
are expected to cancel out in the process of ensemble 
averaging. In order for the two orbits not to be related 
by a symmetry transformation we only need all the K r 
to be nonvanishing. 

Our analysis showed that the behavior of the remaining 
probability density at long times should depend strongly 
on the relative action phase difference between the two 
orbits. This phase difference can be easily controlled by 
adjusting (<j>b — <\> a = NR^/ty- Fixing at a non- 
zero value which produces fa — <j) a = (mod27r), we are 
free to vary K 3 and the three coefficients K' r consistent 
with the constraints 2 - V" > and 2 + T" > (which 
are sufficient to ensure hard chaos). We obtain then, for 
N = 80 and T = 0.1, the results shown in Fig. ||a. Clearly 
the remaining probability is very strongly suppressed on 
the orbit on the left, where the lead is located [as we 
saw in the previous ection, the probability to remain ex- 
actly on the orbit falls off faster with time than prob- 
ability elsewhere, so the numerical value of the relative 
suppression there will be time-dependent]. We also see 
mild density suppression on the two orthogonal invari- 
ant manifolds of this orbit. The phenomenon we want 
to focus on here, though, is the suppression we observe 
on the orbit on the right side of Fig. ||a. The observed 
suppression factor right on the periodic orbit at (1/2,0) 
is 0.37, compared with the predicted value 0.45; again 
the discrepancy may possibly be attributed to finite size 
effects. As expected, we also observe probability suppres- 
sion along the manifolds of this second orbit. 

We now consider the opposite case, where the two or- 
bits are exactly out of phase (fa — 4>a — tt (mod27r)). 
Adjusting K2 appropriately, we again perform ensemble 
averaging over the other parameters and obtain the re- 
sults in Fig. ||b. The same suppression is still seen along 
the orbit containing the lead and its invariant manifolds. 
However, we now see, as expected, an enhancement of 
the remaining probability density near the orbit on the 
right. The enhancement factor on that orbit itself is 1.78; 
the theoretical prediction is 1.67. 

The predicted relative intensity at the orbit (1/2,0) 
(Eq. ^7) given a lead at (0,0) is plotted in Fig. | as a 
function of the action phase difference between the two 
orbits. [One could imagine obtaining such a plot in a 







physical system by tuning a weak magnetic field which 
had little effect on the classical dynamics but did change 
the relative phase between two orbits enclosing different 
amounts of flux. Alternatively, if the periods of the two 
orbits differed, the orbits could be observed to go in and 
out of phase with one another as one changed the energy 
range in which the resonances were populated.] The sup- 
pression and enhancement factors in this case never get 
very far from unity, due to the relatively large value of 
the instability exponent A chosen for our example. Val- 
ues obtained numerically, as described in the preceding 
paragraphs, are given for comparison, along with statisti- 
cal error bars. The dashed line at 1 represents the RMT 
prediction. 

V. CONCLUSION 

We have concentrated throughout on leakage through a 
single decay channel which is very well localized in phase 
space (leak area much less than h). One would like to 
understand more generally possible non-RMT effects for 
multi-channel leads as well as for leads wide enough to 
produce overlapping resonances. A detailed treatment 
of such effects is outside the scope of the present work. 
However, one possible extension turns out to be rela- 
tively straightforward, and the results suggest that the 
most interesting scar effects are already captured in the 
single channel analysis. Specifically, consider M decay 
channels, each with a slow decay rate, as in Eq. |(| If the 
sum of these classical rates is small compared to the level 
spacing, the resulting resonances will be non-overlapping, 
and a perturbative approach to the problem is valid. In 
RMT, the probability to remain at long times is given 
by a product of factors associated with each of the leads 
(Eq. ||). In the case where one of the channels happens 
to be close to a short periodic orbit, the intensity distri- 
bution giving rise to the corresponding factor only will 
be affected, producing the same overall enhancement fac- 
tor < > obtained previously. The total probability 
P rom now falls off much more quickly with time than 
in the single-channel case, reflecting the fact that it is 
much harder to find a resonance which is slowly decaying 
through all of the channels. However, the essential obser- 
vation of an enhanced probability of finding very narrow 
resonances survives. Clearly the analysis is more com- 
plicated for correlated decay channels, and for the case 
of larger total decay rates where the resonances become 
overlapping. Still, periodic orbits by their very nature 
produce quantum effects in phase space regions of size 
h surrounding the orbit; thus our intuition tells us that 
no fundamentally new scar effects are expected in most 
cases for multi-channel leads. 

We have seen that an analysis of short time classical 
motion in a chaotic system can shed much light on quan- 
tum behavior on the scale of the decay time, which is 
much larger than the Heisenberg time and every other 



time scale in the system. This is somewhat counterin- 
tuitive, as the narrow resonance regime is by its very 
nature non-classical and is associated with the very long- 
time behavior of the system. Even though in some cases 
(e.g. in the presence of strong diffraction or caustics) 
semiclassical methods may not be sufficient to predict 
the properties of individual high-energy quantum chaotic 
wavef unctions, they are still very powerful for making 
statistical predictions of the sort described in this paper. 
We note that a small change in the system potential, 
or the presence of a few impurity scatterers may com- 
pletely change the character of individual resonances in 
an open system, making comparison with exact semiclas- 
sical wavefunctions futile. On the other hand, the scar 
theory predictions, which concern statistical properties 
such as the distribution of resonance widths, are robust 
to such changes in the details of the system, as long as 
the dynamics of the first few bounces is known. In a 
situation where the classical dynamics of the quantum 
system under study is not known reliably even for short 
times, one could use methods similar to those described 
here to search for the short unstable periodic orbits. This 
can be done by moving the position of the lead, or, more 
practically in many situations, by adjusting some system 
parameter which changes the classical dynamics of the 
system. For a given classical system and lead position, 
one then sweeps through a weak magnetic field or some 
other parameter not affecting the classical dynamics, and 
searches for a large fraction of very narrow resonances, 
occurring periodically in the magnetic field strength. 

Clearly the ideas described here can also be extended 
to study two-lead systems, where properties such as con- 
ductance peak height distributions can be analyzed. In 
analogy with the present work, the results will strongly 
depend on whether one or both of the leads is located on a 
short unstable classical orbit. Where the leads are found 
on two different short periodic orbits, the phase difference 
between them can be varied to produce an enhancement 
or suppression of the average conductance, as suggested 
by the results of the present work. Even stronger effects 
can be observed if the leads are located on the same peri- 
odic orbit, or if the two orbits are related to each other by 
a symmetry of the system. These issues and other related 
questions are addressed fully in a forthcoming paper [Q . 
Certainly much more work needs to be done generally in 
understanding classical dynamics effects on the quantum 
properties of open chaotic systems. 
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FIG. 1. Smoothed local densities of states Sn n (E) are plot- 
ted as a function of energy on a periodic orbit of instability 
exponent A = 0.20 (solid curve) and on an orbit with A = 0.15 
(dashed curve). The mean resonance width for a lead placed 
on such a periodic orbit will be proportional to Snn(E). We 
observe the peak at the EBK quantization energy E = 
(Eq. |l6| ) which scales as A -1 , the exponential decay between 
£ = and £ = 7T (Eq. and the minimum at the anti-EBK 
energy E = ir, which is exponentially small in A (Eq. |l^). The 
RMT prediction Sn n (E) = 1, which is applicable away from 
any short periodic orbit, is plotted as a dotted line. 



FIG. 2. The probability to remain in the open quantum 
system is plotted as a function of scaled time t' = T c \t. The 
classical prediction exp(-t') is shown as a dotted curve. The 
quantum probability to remain for a generic lead location 
(plusses) compares well with the RMT prediction 1/(1 + t') 
(dashed curve). For a lead placed on a short periodic orbit 
with instability exponent A = 0.96, we obtain the enhanced 
long-time probability to remain (squares), which agrees with 
the scar theory prediction < S^ > /£' (solid line). The sys- 
tem size used for obtaining the data is N = 120 and the decay 
rate per step in the exit channel is T — 0.1. 



FIG. 3. The remaining probability density after time 
t' = 2 • 10~ 3 , as a function of the position of the lead. In 
(a), the numerical data is presented for an ensemble of sys- 
tems of size N = 30; in (b) we show the theoretical prediction 
< > ft' . At the center of each plot is an unstable pe- 
riodic orbit of exponent A = 0.96: for a lead placed at that 
position (white spot) the probability at long times is enhanced 
by a factor of 11 over the same probability for a generic lead 
position (black background). Notice the anisotropy: more 
enhancement at long times is predicted (and observed) when 
the displacement of the lead away from the periodic orbit is 
along one of the invariant manifolds. 
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FIG. 4. The long-time enhancement factor of the probabil- 
ity to remain in a system when the lead is placed on a periodic 
orbit is plotted as a function of the instability exponent of 
the orbit. Data is shown for N = 120 (plusses) and N = 240 
(squares). The N — > oo theoretical prediction < 5,7 > is 
shown as a solid curve. We see the exponential increase in 
the probability to remain as the exponent A decreases (the 
A — > asymptotic form is given in Eq. ^l|). For large A, the 
enhancement factor converges to 1, the RMT prediction. 

FIG. 5. The remaining probability distribution at very long 
times is shown for an ensemble of systems (of size N = 80) 
all having two short periodic orbits in common, each with 
instability exponent A = 1.76. In both cases, the lead is 
centered on the periodic orbit on the left side of the plot. The 
remaining probability is strongly suppressed on that orbit and 
less so on its invariant manifolds. The action of the orbit on 
the right is chosen to be in phase with the first one in case 
(a), so that probability there is also suppressed, and exactly 
out of phase with it in case (b), leading to an enhancement of 
the probability density on the second orbit and its invariant 
manifolds. See next figure for quantitative comparison with 
the theory. 

FIG. 6. The predicted relative intensity at long times 
on a periodic orbit other than the one containing the lead is 
plotted as a function of the relative phase between the two 
orbits (solid curve). For reference, the phase-space averaged 
intensity is plotted as a dashed line. Intensity suppression 
is predicted and observed when the orbits are in phase, and 
enhancement is seen when the orbits are exactly out of phase 
[compare (a) and (b) in previous figure] . The data is collected 
for an ensemble of systems with N = 80, and each orbit has 
instability exponent A = 1.76. 
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